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CONNECTING SUBSTRUCTURE IN GALAXY CLUSTER CORES AT Z=0.2 WITH CLUSTER ASSEMBLY HISTORIES 
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ABSTRACT 

We use semi-analytic models of structure formation to interpret gravitational lensing measurements of sub- 
structure in galaxy cluster cores (/?<250/z _1 kpc) at z = 0.2. The dynamic range of the lensing-based substruc- 
ture fraction measurements is well matched to the theoretical predictions, both spanning /sub^ 0.05-0.65. The 
structure formation model predicts that / su b is correlated with cluster assembly history. We use simple fitting 
formulae to parameterize the predicted correlations: A90 = T9o+a9olog(/ su b) and A50 = T5o+ct!5olog(/ SU b), where 
A90 and A50 are the predicted lookback times from z = 0.2 to when each theoretical cluster had acquired 90% 
and 50% respectively of the mass it had at z = 0.2. The best-fit parameter values are: 090 = (-1 .34 ± 0.79) Gyr, 
r 90 = (0.31 ±0.56) Gyr and a 50 = (-2.77 ± 1.66) Gyr, r 50 = (0.99 ± 1.1 8) Gyr. Therefore (i) observed clusters 
with / su b<; 0.1 (e.g. A 383, A 1835) are interpreted, on average, to have formed at z J> 0.8 and to have suffered 
< 10% mass growth since z — 0.4, (ii) observed clusters with / su b^ 0.4 (e.g. A 68, A 773) are interpreted as, on 
average, forming since z — 0.4 and suffering > 10% mass growth in the ~ 500 Myr preceding z = 0.2, i.e. since 
z = 0.25. In summary, observational measurements of / su b can be combined with structure formation models to 
estimate the age and assembly history of observed clusters. The ability to "age-date" approximately clusters in 
this way has numerous applications to the large clusters samples that are becoming available. 

Subject headings: gravitational lensing — cosmology : dark matter — galaxies:clusters:individual (A 68, A 209, 
A267, A383, A 773, A963, A 1763, A 1835, A2218, A2219) 



1. INTRODUCTION 

The mass growth of clusters is sensitive to the dark energy 
equation of state parameter w, the matter density of the uni- 
verse VL m and the normalizat i on of the ma t ter power spectrum 
o-g (e .g. lEvrard et all 119931: ISmith et al.1 120031: lMantz etal.1 
120071) . Clusters are inferred to grow hierarchically via the 
ingestion of smaller dark matter halos (that host galaxies) 
into the more massive parent halo (cluster). The structure of 
galaxy clusters, specifically the internal substructure of clus- 
ters, therefore contains a wealth of cosmological information, 
including possible clu es about the physics o f the dark mat- 
ter particle itself (e.g. iNataraian et al.1 120021) . From an as- 
trophysical point of view, the mass growth of clusters brings 
new (generally gas rich) galaxy populations into clusters (e.g. 
iMoran et al.ll2007l) and may lead to shock-heating of the intr- 
aclust er medium and/ or disruption of cooling in cluster cores 
(e.g. lPoole et al.ll2008l) . Reliable measurement and intepreta- 
tion of cluster substructure is therefore of broad interest. 

The most direct way to detect substructure within clus- 
ters is via gravitational lensing. Group-scale substructures 
within individual clusters were detected in ear ly ground- 
based strong-lensing studies of individual clusters (Pel lo et al.l 
I199U iKneib et al.1 11993L 1 1995b and subsequently measured 
to high precision u si ng Hubble Space T elescope (HST) data 
(Knei b et al.lfl996l) . ISmith et al l (120051 - hereafter Sm05 - 
see © then measured the structure of a sample of 10 clus- 
ters at z~0.2. In this article we use iTaylor & Babulf s (12004! 
- hereafter TB04) semi-analytical models of structure forma- 
tion to interpret Sm05's cluster substructure measurements, 
as a means of exploring lensing-based substructure measure- 
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ments as a quantitative probe of cluster age and assembly his- 
tory. 

We summarize Sm05 and TB04 in ^2] and ^respectively, 
and then synthesize observations and theory in $4] We dis- 
cuss caveats in ^5] and summarize our conclusions and dis- 
cuss future prospects in §|6] We assume Ho=70kms~ 1 Mpc~ I , 
17ivi=0.3, f2A=0.7 and erg = 0.9 throughout. The lookback time 
from z = to z = 0.2 is f, = o.2 = 2.44Gyr in this cosmology. 

2. SUMMARY OF OBSERVATIONAL RESULTS 

Sm05 investigated the projected mass and structure of ten 
X-ray luminous (L x > 4 x 10 44 ergs _1 , 0.1 -2.4keV) cluster 
cores at 0.17 < z < 0.25 (Table [TJ. Hubble Space Telescope 
( HST) AVFPC2 imag ing and ground-based spectroscopy of 
gravitational arcs (ISmith et alfcoOll l2002t ISand et al.ll200l 
Richa rd et al.1 12007F were used to characterize the strong 
and weak gravitational lensing signal of each cluster core. 
The lensing signals were then used to constrain a detailed 
parametrized model of the projected mass distribu tion in each 
cluster core fol lowing Kneib (19931- see also Kne ib et all 
1996 and Smith 2002J). Each lens model includes mass com- 
ponents that account for both the underlying dark matter dis- 
tribution in the cluster (cluster/group-scale mass components) 
and the cluster galaxies down to Lk>0.1L* k . For clarity, we 
refer to the main central cluster dark matter halo as the cluster- 
scale mass component, and all other massive substructures as- 
sociated with infallen clusters and/or groups as group-scale 
mass components. The typical number of galaxy-scale mass 
components in each lens model was 30. 

The HST data probe out to a typical projected cluster- 
centric radius of R = 250 /r'kpc at the cluster redshifts. Sm05 
measured the cluster substructure within this radius. Here 
we define / su b as the fraction of mass associated with sub- 
structures, / sub =M sub /M tot = 1 -M cen /M to t, where M cen /M tot 
is the central mass fraction used by Sm05. For this purpose, 
"substructures" include group-scale masses and galaxy-scale 
masses, with the exception of the brightest cluster galaxy 
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TABLE 1 

Observational Substructure Measurements 



Cluster 


Redshift 


Mux 


Ndm 


/sub 


Strong Lensing? 


Sm05 Classification 


A 50 


Ago 






(10 14 M Q ) 










(Gyr) 


(Gyr) 


A 383 


0.188 


2.6±0.1 




0.06±0.01 


Confirmed 


Undisturbed uni-modal 


4.3±2.4 


2.0±1.1 


A 963 


0.206 


2.4 ±0.2 




0.06±0.01 


Confirmed 


Undisturbed uni-modal 


4.3±2.4 


2.0±1.1 


A 1835 


0.253 


4.1 ± 1.1 




0.06±0.01 


Unconfirmed 


Undisturbed uni-modal 


4.3±2.4 


2.0±1.1 


A 267 


0.230 


1.9 ±0.4 




0.08±0.01 


Unconfirmed 


Disturbed uni-modal 


4.0±2.2 


1.8±1.0 


A 1763 


0.288 


1.5±0.8 




0.15±0.05 


No? 


Disturbed uni-modal 


3.2±1.8 


1.4±0.9 


A 209 


0.209 


1.1 ±0.5 




0.18±0.06 


No? 


Disturbed uni-modal 


3.0±1.7 


1.3±0.8 


A 2219 


0.228 


2.4±0.1 


3 


0.20±0.01 


Confirmed 


Disturbed multi-modal 


2.9±1.7 


1.3±0.8 


A2218 


0.171 


4.0 ±0.1 


2 


0.27±0.01 


Confirmed 


Disturbed multi-modal 


2.5±1.5 


1.1±0.7 


A 68 


0.255 


3.1 ±0.1 


2 


0.36±0.01 


Confirmed 


Disturbed multi-modal 


2.2±1.4 


0.9±0.7 


A 773 


0.217 


3.6±1.2 


3 


0.61±0.20 


Unconfirmed 


Disturbed multi-modal 


1.6±1.3 


0.6±0.6 



(BCG). This is because BCG mass is degenerate with cluster- 
scale mass in the lens models, and, in any event, Sm05 found 
no evidence for BCGs being offset from the center of the 
cluster-scale mass components. The clusters are listed in 
order of increasing / su b in Table Q] together with Ndm, the 
number of cluster/group-scale mass components in each lens 
model. As expected, A^dm and / su b are correlated, however 
there is a factor of 3 spread in / su b for sub-samples of clus- 
ters with Num= 1 and /Vdm> 1- Note, that to aide compari- 
son of Sm05's lens modeling results with predictions, the val- 
ues of / su b in Table Q] have been adjusted upwards to take ac- 
count of the fact that the "L"-shaped observed WFPC2 field 
of view (~ 5 arcmin 2 ) covers just ~ 50% of a circle of radius 
250/r'kpc (~ lOarcmin 2 ). The WFPC2 observations were 
originally designed to include the likely group-scale substruc- 
tures in the cluster cores. The adjustments to / su b therefore 
account statistically for the galaxy-scale masses not included 
in Sm05's analysis, leaving Ndm unchanged. 

Cluster-to-cluster differences in / su b arise for two reasons: 
(i) group-scale mass components in the lens models that are 
associated with massive, likely in-falling structures such as 
groups of galaxies, and (ii) cluster galaxies that are associated 
both with the central cluster dark matter halo (and are pre- 
sumably virialized) and with the in-falling structures. Sm05 
broadly interpreted measurements of / su b as indicating the 
merger history of clusters, however quantitative conclusions 
on cluster assembly and age were impossible without theoret- 
ical models. 

3. SUMMARY OF THEORETICAL MODELS 

TB04's semi-analytic model of halo evolution provides a 
fast way of generating a large number of model halos for 
comparison with observations, and agrees well in its predi- 
tions with self-consistentn-body simulations (Tay lor & Babull 
2005). The model includes two main components, a merger- 
tree code for determining the assembly history of an individ- 
ual cluster, and an analytic description of how the main halo 
and the merging subcomponent evolve after each merger. We 
will summarize these components here, and refer the reader 
to TB04 for full details. 

The merger tree describing the assembly of a sin- 
gle cluster is generated randomly using the algorithm of 
ISomerville & Kolattl (1 19991) . Starting from a halo of speci- 
fied mass at z = 0, the algorithm chooses an interval to step 
back in redshift and picks progenitors at this r edshift follow- 
ing extended Press-Schechter merger statistics dLacev & Cold 



119931) . Iterating produces a complete history of the mergers 
through which the final object assembled, down to some lim- 
iting mass resolution and back to some redshift. In this article 
we use a set of 1000 merger trees whose final masses at z = 
are randomly drawn from the massive end (> 5 x 10 14 M Q ) 
of a halo mass function. For the range of concentration 
parameters derived for our halos (as explained below), this 
produces projected masses within a cluster-centric radius of 
R = 250 /r'kpc (integrating out to the virial radius along the 
line of sight) in the range 0.4-8 x 1O 14 M , similar to those 
in the observed sample (Table [TJ. The merger trees for these 
systems are followed back to a redshift of 50, or until they 
drop below the resolution limit of lO^x the final mass. 

Given a merger tree, the model of TB04 selects the most 
massive progenitor in the most recent merger and traces its 
history back, selecting the most massive progenitor at each 
time step. This object is considered the "main" system, and 
is modelled as a spherical halo with a radial density pro- 
file similar to the " universal" profile found in simulations 
dNavarro et al.| [l997). To account for some of the centrally 
concentrated (baryonic) m ass in the cl u ster, w e use the fit- 
ting formula proposed by iMoore et al.l d 1998b . which has a 
steeper r" 1 5 central cusp. In an update to the model presented 
in TB04, in this work we use the whole mass assembly his- 
tory of the system, rather than just its instantaneous mass, to 
derive the value of its concentration parameter at each redshift 
step, following Wec hsler et al.1 (120021) . Each merging subhalo 
is added to this main system as a smaller spherical object with 
a Moore profile, and its subsequent orbit and mass-loss his- 
tory are determined using the analytic descri ption of dynami- 
cal fric tion, tidal stripping and tidal heating in Tay lor & Babull 
(1200 ll) . Subhalos are tracked until they are disrupted either 
by repeated tidal stripping or by passing within 0.01 ■/?,,,> of 
the centre of the main potential. The model also includes sim- 
ple treatments of sub-substructure (subhalos within subhalos), 
correlated orbits in infalling groups of objects, and collisions 
and encounters between objects. See TB04 for a detailed ex- 
planation of these elements. 

For each merger tree, we use the semi-analytic model to de- 
termine how much substructure exists at z = 0.2, and calculate 
/ su b as the fraction of the mass projected within 250 /r'kpc of 
the halo centre that is contained in the 30 most massive sub- 
halos in the system (to match the observations - see We 
exclude from the final analysis any substructure which ends 
up within 5% of /?„> of the halo centre, since in a real cluster 
most such objects would interact strongly and merge rapidly 
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FIG. 1 . — Predicted distribution of substructure fraction versus (90 (left) and 
?50 (right). Note that the lookback time from the z. = to z. = 0.2 is 2.4 Gyr 
(dashed line). The grey open points show clusters from TB04's model - they 
span a similar dynamic range to the observed clusters (see Table[T}. The big 
black points show the mean prediction in equally spaced logarithmic bins; 
the error bars show the lc scatter in each bin. The solid black line shows the 
best fit relations described in the text (@ 

with the dominant central galaxy. We also record several in- 
dicators of the mass assembly history of the main system in 
each merger tree, including £90 and tso, the lookback times by 
which it had assembled 90% and 50% respectively of the mass 
it has at redshift 0.2. These thresholds are chosen to probe the 
recent infall history and cluster age respectively. We show the 
resulting / S ub - tgo and / su b - ?5o distributions in Fig. Q] The 
observed clusters span a dynamic range of 0.06</ su b<0.62 
(Table [T), which is comparable with that of the TB04 theo- 
retical models, and consistent with other theoretical estimates 
of cl uster substructure ("e.g. [Diemand et al.ll2004t Ipao et alJ 
12004 iNagai & Kravtsovll2005l:lvan den Bosch et al.ll2005l) . 

4. RESULTS 

Fig. Q] reveals that the TB04 model predicts that / su b and 
tgo are correlated, as are / S ub and t$o. However, at / su b< 0.2 
the distributions of tgo and fg are bi-modal - some theoretical 
clusters exhibiting both a low substructure fraction and recent 
significant infall. The unphysical location of these clusters in 
Fig. 1 is due to a timing difference between the axes - tgo is 
sensitive to mass growth within the virialized region of the 
cluster (approximately a sphere of diameter ~ 3Mpc); / su b is 
sensitive to substructures in a cylinder of diameter 500 /i"'kpc 
through the center of the cluster. If an infalling structure 
has crossed the virial radius but not yet entered the cylinder, 
then it may appear in the "spike"of clusters at / su b< 0.2 and 
tgo< 3Gyr (or ?5o< 3Gyr). We exclude these clusters from the 
calculations described below. 

We have therefore confirmed and considerably extended 
Sm05's qualitative interpretation of the lensing results: clus- 
ters with high substructure fractions have (i) suffered more 
pronounced recent infall than clusters with low substruc- 
ture fractions, and (ii) formed at later times than clusters 
with lower substructure fractions. To quantify this, we 
fit the following simple formulae to the theoretical data: 
A 90 =tg -t z= o.i = r 90 + Q!9olog(/ su b) and A 50 =f 5 o -t z= o.2 = t 50 + 
q?50 log(/ su b), where ago and ago parameterize the dependence 
of A90 and A50 respectively on / su b; Tgo and Tgo are the in- 
tercepts at / su b= 1. We first bin up the individual theoretical 
data points in bins of width Alog(/ S ub) = 0.2, and then weight 
each bin in the fit by the reciprocal of the sample variance in 
the respective bin. The best fit parameter values obtained in 
this way are: ago = (-1 .34 ± 0.79) Gyr, r 90 = (0.3 1 ± 0.56) Gyr 
and a 50 = (-2.77 ± 1.66) Gyr, r 50 = (0.99 ± 1.18) Gyr. The 



Tgo and T50 parameters are measurements of the timing dif- 
ference between the axes in Fig. Q] discussed above, and are 
well-matched to the infall timescale of ~ 0.5 - 1 Gyr. 

The best-fit models were then used to interpret quantita- 
tively the observed substructure fractions from Sm05 listed in 
Table Q]- i.e. to estimate the age (A50) and recent infall his- 
tory (A90) of each cluster - see TableQ] The uncertainties on 
A90 and A50 quoted in Table Q] incorporate errors on / su b and 
on the best-fit parameter values derived above, and are dom- 
inated by the scatter on / su b - £50 and / su b - £90- The typical 
uncertainty on cluster age is ~ 1.7 Gyr, with ages spanning 
~ 1 -4 Gyr . Cluster with the largest substructure fractions, 
i.e. /subjj 0.4 formed within the ~ 2Gyr preceding z = 0.2, i.e. 
since z <J 0.4, and had not assembled 90% of the mass they 
had at z = 0.2 until z ~ 0.25, i.e. just ~ 500Myr before z = 0.2. 
In contrast, clusters with the lowest subtructure fractions, i.e. 
/sut>< 0.1 formed > 4Gyr before z = 0.2, i.e. at z ;> 0.8. and 
then went on to assemble > 90% of their z = 0.2 mass within 
the ensuing ~ 2 Gyr, and suffered negligible mass growth in 
the ~ 2 Gyr prior to z = 0.2, i.e. since z — 0.4. Therefore we 
interpret clusters with the lowest / su b as being, on average, 
almost fully (~ 90%) assembled by the time that the clusters 
with highest / su b had barely formed (i.e. assembled half of 
their mass). 

5. CAVEATS 

TB04 make some simplifying assumptions that may affect 
our results: (i) dark matter halos are assumed to be spherical, 
whereas real dark matter halos are likely triaxial; (ii) all mat- 
ter is treated as collisonless, thus ignoring baryonic physics 
including adiabatic contraction due to gas cooling and the 
dynamical effects of galaxies; (iii) matter outside the clus- 
ter virial radius is ignored; in contrast, lensing is sensitive to 
all the matter along the line of sight, including chance fore- 
ground and background projections. We will consider these 
complications in more detail in future work. 

The main systematic uncertainty on the observed substruc- 
ture fractions is the completeness of Sm05's lens models as a 
function of sub-halo mass. A detailed study of this issue will 
be published in the future (Hamilton-Morris et al., in prep.). 
Here we identify how many of Sm05's models may suffer 
from incompleteness. Incompleteness is most likely to arise 
from features in the dark matter distribution not identified in 
the luminous properties of the clusters within R < 250 /r'kpc, 
including group-scale dark matter halos in which the cluster 
galaxies may be embedded. Strong lensing has been detected 
in eight of the ten observed clusters (Table [TJ; group-scale 
dark matter halos are detected in three of the eight by the mea- 
surable way they alter the appearance of strongly-lensed arcs 
compared to simpler mass distributions. Clusters in which 
strong-lensing has not been detected are therefore of greatest 
concern, as pointed out i n Sm05's discus si on of t he filament 
feeding A 1763; see also Mer curio et alj s ([2003) discussion 
of the complex dynamical structure of A 209. Gross substruc- 
ture incompleteness may therefore be a concern for ~ 20% of 
the observed sample. 

We also consider whether there is any evidence in the re- 
spective / su b distributions for systematic differences between 
Sm05's X-ray selected sample and TB04's mass-selected 
sample. To address this we compare the cumulative and dif- 
ferential / su b distributions in Fig. [2] A KS test obtains a prob- 
ability of 20% that the observed and theoretical distributions 
are drawn from the same parent distribution, i.e. formally ev- 
idence of systematic differences between the samples at 1.3cr 
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FIG. 2. — Distribution of substructure fractions from the observed sample 
(black solid line) and from the theoretical model (red dashed). The observed 
clusters show a possible excess of very low and very high substructure sys- 
tems - this is discussed in iJ5] 

significance. The small observed sample size of just 10 clus- 
ters is clearly a limiting factor here. The right panel of Fig. [2] 
helps somewhat, as it highlights more clearly the possible bias 
of the observed sample to very high and very low values of 
/ Sll b. Such a bias would be plausible because cluster-cluster 
mergers and cool cores both suffer an excess of X-ray flux at 
a fixed mass over a non-merging, non-cool core cluster. A 
much larger observational sample is required to investigate 
this issue further. 

6. CONCLUSIONS AND FUTURE PROSPECTS 

We have combined theoretical models of structure forma- 
tion (TB04) with gravitational lens models of galaxy clusters 
(Sm05) to explore how measurements of cluster substructure 
from lensing observations can be interpreted in the context 
of the age and assembly history of clusters. The main re- 
sult is that / su b, the fraction of cluster mass within a projected 
cluster-centric radius of R = 250 /i _1 kpc associated with sub- 
structure (galaxies and group-scale halos), as can be mea- 
sured from lensing data, is predicted to be strongly corre- 
lated with the age and recent mass growth of galaxy clus- 



ters. We fitted the the following simple formulae to the the- 
oretical data to quantify the predicted behavior in a conve- 
nient form: A 90 =f 90 ~t z =o.2 = t 90 + a 90 log(/ su b) and A 50 =f 90 
-t z =o.2 = T50 + CH50 log(/ su bX where t% and tso are the look- 
back times at which a cluster had acquired 90% and 50% 
of its mass at z = 0.2. The best fit parameter values are: 
a 90 = (-1.34±0.79)Gyr, r 90 = (0.31 ± 0.56)Gyr and a 5 o = 
(-2.77 ± 1.66) Gyr, r 50 = (0.99 ± 1.18) Gyr. Low-/ sub clusters 
Cfsub<; 0.1; e.g. A 383, A 1835) are therefore interpreted as, on 
average, having formed at z<;0.8 and having suffered <10% 
mass growth in the 2 Gyr preceding z = 0.2, i.e. since z — 0.4. 
In contrast, high-/ su b clusters (/ su b^ 0.4; e.g. A 68, A 773) 
are interpreted, on average, to have formed just ~ 2 Gyr be- 
fore z = 0.2, i.e. since z — 0.4, and suffered 10% mass growth 
in the ~ 0.5 Gyr preceding z = 0.2, i.e. since z~0.25. 

Our synthesis therefore demonstrate that lensing-based 
measurements of / su b can be combined with semi-analytic 
structure formation models to estimate the average age and 
assembly history of observed clusters. This suggests numer- 
ous avenues for further exploration, including: (i) expan- 
sion of the observed samples by at least an order of magni- 
tude, (ii) calibration of the completeness of the lensing-based 
mass function of sub-halos in clusters, (iii) investigation of 
how lensing-based cluster age and assembly history estimates 
might allow new cosmological constraints to be derived, for 
example, on the dark energy equation of state parameter w, 
(iv) analysis of cluster galaxy populations and cluster scaling 
relations as a function of cluster age. 
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